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Abstract 

We discuss the possibility of formation and subsequent detection of a supersymmetric bound 
state composed of a slepton-antislepton pair at the next linear collider. The Green function 
method is used within a non-relativistic approximation to estimate the threshold production cross- 
section of the 2P bound state. The parameter space of Gauge Mediated Symmetry Breaking 
(GMSB) models allow a particular scenario in which a charged slepton (e/j, JIr or ri) is the NLSP. 
Within this scenario the produced 2P bound-state decays, through a dipole transition, into the 
IS* ground-state with branching ratio ~ 100% emitting a very soft (~ 1 MeV) photon which goes 
undetected. The spectroscopy of the IS'-state shows that it decays into two photons with Br « 0.5 
up to niNLSP ~ 1 TeV. Thus NLSP sleptonium threshold production gives rise to the signal 
e^e“ —> 2P —> 15 -|- ^^soft 7 ” —> 77 which when compared with the standard model two-photon 
process (e^e“ ^ 77 ) has a statistical significance (55=signal/noise) which, at an energy offset 
from threshold of F = 20 GeV, goes from SS = 11 to SS = 2 when the mass of the NLSP ranges 
in the interval [100, 200] GeV. 

PACS numbers: 12.60.Jv, 11.10.St, 14.80.Ly 
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I. INTRODUCTION 


Despite the enormous success of the standard model (SM) of particle interactions based 
on the gauge group SU{3)c ® SU{2)l ® U{1)y a large portion of current research efforts in 
the held of fundamental interactions is devoted to the study of signatures of physics beyond 
the SM. Between the possible alternatives to the standard model its minimal supersymmetric 
extension (MSSM) is one of the theories which has been extensively studied from both the 
purely theoretical and phenomenological aspects [1]. 

Supersymmetry, is the symmetry which relates fermions to bosons and, while appealing 
because of its potential of solving the hierarchy problem of the standard model, must of 
course be broken, and one of the major issues in supersymmetric theories is the pattern 
of supersymmetry breaking. There exist various possibilities to break supersymmetry. In 
mSUGRA models it is assumed that supersymmetry is broken at the Planck scale Mp and 
is transmitted to the low-energy sector by gravitational interactions only [2]. In gauge 
mediated symmetry breaking (GMSB) models supersymmetry is broken at relatively lower 
energy scales and is mediated by gauge interactions [3]. Other possibilities consist of anomaly 
mediated symmetry breaking (AMSB). Here the SUSY-breaking occurs also in a hidden 
sector but it is transmitted to the visible sector via the super-Weyl anomaly [4]. 

Of course at this stage of experimental and theoretical investigations it is important to 
study in detail the phenomenological consequence of each scenario. Object of this work is 
the GMSB scenario whose theoretical basis and phenomenology has already been described 
in the literature, (see [3] and references therein). In particular we recall that GMSB models 
are characterized by a rather peculiar superparticle mass spectrum. Indeed in these models 
the almost massless goldstino/gravitino, G, is the lightest supersymmetric particle (LSP) 
and in addition to the neutralino Xi next-to-lightest supersymmetric particle (NLSP) case, 
large fractions of the parameter space offer the possibility of having a charged NLSP which 
can be either the ri or ip {I = e, /i) . 

In [5] the phenomenology of the production of a pair of NLSP, within GMSB models, 
was discussed in the context of the Gern LEP2 collider. In particular it was found that in 
the stau NLSP or charged slepton NLSP cases typical signatures of NLSP pair production 
are If, or i~^i~]f where the missing energy. If, is due to the decay of the NLSP to the 

almost massless gravitino, G, (LSP). 
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The object of this study is to propose a novel type of signature for the GMSB slepton 
NLSP scenario, assuming that the next LC will operate (at least initially) at energies approx- 
imatively in the neighborhood of the threshold of NLSP pairs. In this case if or rj are the 
NLSP then the threshold production of sleptonium, a bound state of a slepton-antislepton 
pair, (smuonium or stauonium) in e’''e“ collisions should be considered. It turns out that 
such bound state can give rise to the interesting signature of two photons and practically 
no missing energy. This signature has not been considered previously. Indeed the produced 
2P state decays with branching ratio ~ 100% to the IS* state and a photon whose energy 
is related to the difference of the energy levels E 2 P — Eis ~ 1 MeV. In turn the IS* state 
decays mainly to two photons: 


e’''e — 2P —IS* + soft 7 

U 77 ( 1 ) 

The threshold production cross section of the 2P state in e’''e“ collisions is computed using 
the Green function method within a non relativistic approximation. The decay widths of all 
open channels of the 2P and IS* states are also provided. It is concluded that the two-photon 
signature, when compared with the SM process e’''e“ —> 77 , has a statistical signihcance 
{SS = Signal/Noise) which ranges in the interval from 11 — 3 when the mass of the slepton 
{piNLSp) is between 100 and 300 GeV. 

The plan of the paper is as follows: section II provides a review of the GMSB supersym¬ 
metric model; in section III is discussed the criterion for formation of the supersymmetric 
bound state (sleptonium); section IV presents details of the decay channels of the 2P and 
IS* bound states; section V contains a description of the Green function method to estimate 
the threshold cross-section for the production of the 2P bound state; section VI presents 
a discussion of the statistical signihcance of the two-photon signature; hnally section VII 
presents the conclusions. 


II. GMSB MODELS 

GMSB models of symmetry breaking are perhaps the most promising alternative to the 
SUGRA scheme where SUSY breaking takes place at the Planck mass and is then com¬ 
municated to the low energy sector by gravitational interactions. In GMSB models super- 
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symmetry breaking occurs at relatively lower energy scales and it is mediated by gauge 
interactions [3, 6 ]. One nice feature of these models is the automatic suppression of the 
SUSY contribution to FCNC and CP-violating processes. 

In the simplest version such models are characterized by the introduction of messenger 
chiral superheld which contain quarks 'ipq,'>pq, leptons 'ipij'ipi, scalar quarks q,q and scalar 
leptons i,i (messenger fields). All these particles must acquire very large masses as they 
have not been discovered. They do so by coupling to a gauge singlet chiral supermultiplet 
X via the superpotential: 

W^ess = X2Xei+X3Xqq ( 2 ) 


Then one assumes that the scalar component of X and its auxiliary F-term acquire a vacu¬ 
um expectation value (VEV) respectively denoted (S) and (Fs). Assuming for simplicity 
degeneracy of the couplings (A = A 2 = A 3 ) and absorbing the coupling A into (S') and a/(F 5 ) 
by dehning M = |A(S')| and F = |A(F 5 )| the amount of SUSY breaking in the messenger 
sector i.e. the mass splitting of the scalar messenger states is found to be parametrized as: 


77 ?^ = 

''^mess. fermions ^ 


"^mess. scalars = M ±F 


6m 


2 _ 9 p 

mess, scalars 


(^TyT-mess. scalars ^ 


(if F/M^ < 1) (3) 


SUSY breaking is thus apparent in the messenger sector, and is in turn communicated to 
the low energy sector (MSSM sparticles) through radiative quantum corrections. Gauginos 
(and scalar partners) acquire their masses, at the messenger scale M, through one loop (and 
two loop) Feynman diagrams where virtual messenger particles are exchanged [7-9]: 


Mg{M) = -^Agix) G = U{1),SU{2),SU{3) (4) 

m 

ml{M) = 2A‘^ G^{s)N^{m) (5) 

G,m 

In Eqs. (4,5) x = F/M"^, m labels the messengers and s the MSSM scalar; the functions 
f{x) and g{x), which reduce to ~ 1 when x ^ 0 are explicitly given in [5, 8 , 9]; 
is the Dynkin index of the gauge representation under which the messenger superhelds 
transform, dehned by TrT“T^ = {N^/2)6°'^ the T“ being the generators of the gauge group 
in the representation R^; are the quadratic Casimir invariant of the same gauge group 


1 = (6/5)r2, where Y = Qem - T 3 . 
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representation for the MSSM scalar field in question and, for the N of SU{N)), is defined 
by: = (iV^ — l)/2iVJ and turn out to be simple algebraic 

functions of the gauge couplings and the number and riq of messenger helds, see [5, 9] for 
further details. From Eqs. (4,5) one also deduces that the scale of SUSY breaking felt in the 
messenger sector A = F/M must be in the range 

10 TeV < A < 100 TeV 


in order to have sparticle masses in the range of 100 GeV - 1 TeV. 

A distinctive feature of GMSB models is the fact that the gravitino G may be very light. 
Indeed mg is given by: 


rriQ = 1713/2 


fq f y 

V3M'p ^ VlOOTeVy 


2.4 eV 


( 6 ) 


where Mp = = 2.4 x 10^® GeV is the reduced Planck mass and is the 

fundamental scale of super-symmetry breaking (SSB) which does not coincide with F, the 
scale of SUSY breaking felt by the messenger sector. The ratio F/Fq depends on how 
SUSY breaking is communicated to the messenger sector. If the communication takes place 
via a direct interaction then F/Fq is given by the corresponding coupling constant which 
by imposing perturbativity arguments can be shown to be smaller than 1 [5], thus giving 
Fo > F. If the communication of SUSY breaking takes place radiatively then F/Fq is given 
by some loop factor and thus F/Fq 1. It can easily be shown [5, 10] that Fq is only 
subject to a lower bound which is: 

> A (7) 


which is typically of the order of 10 — 100 TeV. 

Therefore the gravitino G turns out to be the lightest super-symmetric particle (LSP). In 
R conserving supersymmetry all sparticles eventually decay to the gravitino and in order to 
compute the decay widths one needs the interaction Lagrangian in the gravitino held which 
can be computed in the limit of global supersymmetry ( if Mp) as the dominant 

gravitino interactions come from its spin 1/2 component (the goldstino). It is therefore 
a good approximation to describe the gravitino LSP in terms of its spin 1/2 goldstino 


2 (jUy{1) = (3/5)r2 
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component. Goldstino interactions contain derivative couplings suppressed by I/Fq: 

£in. = -jr (MYd.ip - d,G + h.c. (8) 

Sparticle phenomenology (production and decay) is strongly affected by the type of the next 
to lightest super-symmetric particle (NLSP). All sparticles will decay to a cascade leading to 
the NLSP which will in turn only decay to the gravitino G via I/Pq interactions. Depending 
on the values of the parameters, the NLSP can be either the neutralino Xu ^^e stau ri, or 
in restricted regions of the parameter space the sneutrino (P). 

Of particular interest in this work is the case of a charged NLSP. In GMSB models the Ti is 
always the NLSP but in some circumstances it may happen that the mass of the be closer 
to niri than the mass of the tau (m,-). If this is the case the Ir = (cr, JIr) act effectively as 
a NLSP since the decays Ir are kinematically forbidden. This situation is referred 

to as slepton co-NLSP scenario, and it is more precisely dehned by the condition; 

< Min [m^o, -|- nir] ■ (9) 


Within this scenario the sleptonium bound state of a pair of ri or Ir would be the lightest 
SUSY state to he produced in a laboratory. It is therefore interesting to explore throughly all 
possibilities to detect such a bound state at the next linear collider (NLG). In the following 
section we describe the spectroscopy in detail. 

Within this scenario the NLSP i = {cr, jiR, Ti) total width is easily determined from Eq. (8): 


Vl = V{l^^G) 


167rF2 


( y 

VlOOGeV/ 


/100TeV\ 


4 

2 X 10"^ eV 


( 10 ) 


which will be used to establish a criterion for the formation of the bound state. 

It should be noted that Eqs.(4,5) are to be considered as boundary conditions at the 
(high) messenger mass scale M. Low energy values of the parameters are to be obtained 
by running renormalization group equations (RGE) down to the electroweak scale. This 
process must of course ensure proper breaking of the electroweak (EW) symmetry. All 
this is achieved by the public domain code Suspect [9], a software which allows for the 
possibility of choosing between mSUGRA, GMSB and AMSB models in addition to the 
general unconstrained MSSM. With Suspect is possible to perform a scan of the parameter 
space in order to select the scenario that one is interested in. 
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In its simplest version the GMSB model is characterized by a relatively small number of 
parameters: 

F 

M, ^ tan(/3), sign{^i) (11) 

In Fig. 1 we show the result of such a (numerical) study of the parameter space having 
hxed tan/? = 4 and ni = riq = 5. The scan has been performed on the plane M, A. In 
the hgure we show the masses of ^NLSp) and the higgs mass states entering into the 
spectroscopy of the bound state decays to be discussed in sec. IV. The mass of the neutralino 
is shown for completeness to conhrm that it is heavier than 


III. FORMATION SCENARIO 

In this section we will review the creation of the bound state. For the SUSY case, our 
assumption will be that the creation of the bound state does not differ from the Standard 
Model case, as the relevant interaction is again driven by QED and is regulated by the mass 
of the constituent superparticles. A criterion for the formation of bound states we shall 
adopt is that [11, 12] the formation can occur only if the level splitting, which depends upon 
the strength of the interaction among the (s)particles, is larger than the natural width of the 
would-be bound state. It means that the bound state is formed if the following condition is 
satished 

AE 2 P- 1 S > F (12) 

where AE 2 P- 1 S = E 2 P — EiSi and F is the width of the would-be sleptonium, which is twice 
the width of the single slepton F = 2F^, as each slepton could decay in a fashion independent 
from the other. F is not the total decay width of the sleptonium bound state, as it includes 
only the single smuon decay modes and not the annihilation modes. It represents the minimal 
energy level spread necessary for bound state formation. If created, the bound state will in 
turn also have its own annihilation decay modes (as discussed in sec. IV). 

For the case of a scalar bound state {sleptonium), we should consider the Coulombic two- 
body interaction 

U(r) = (13) 

r 

where the coupling a is the usual fine structure constant of QED. With this position we are 
able to compute analytically the energy levels and the wavefunctions within a non relativistic 
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approach. From 


(14) 




mi o? 

4 V? 


one infers that 


AF^ 2 P-is' = —~ 1 MeV 


mi 


16^' LlOOGeVJ 

Contrary to the QCD interaction case [13] the running of the coupling constant value is not 

very important, as the relevant scale given by the Bohr radius, 2/{mia)^ is of 0(1) GeV“^. 


The AE value is thus determined only by the mass of the slepton. This has to be compared 
to the width of the would-be sleptonium T = 2r^. Thus the requirement of formation, is 
obtained inserting Eq. 15 and Eq. 10 in Eq. 12: 


Bound State formation 



( 16 ) 


We must emphasize that formation criterion adopted here for the sleptonium bound 
state is slightly less stringent than the one based on the revolution time, for which no bound 
states exist, if the revolution time, tn = 27rrjv, is larger than the lifetime of the rotating 
constituents, r = l/T, as shown in [12], 

In order to estimate the revolution time, we use the consequences of the virial theorem, 
which reads (T) = —{V)/2 for the average of kinetic and potential energies respectively. 
From the expression for the energy levels Eq. 14 we obtain the average speed of the con¬ 
stituent slepton, = a‘^/{n‘^), that is 


(n) = - (17) 

n 

while the average distance of the constituent is given by 

(r) = 

for a Coulombic potential as in Eq. 13. Combining Eqs. 17 and 18 we compute the revolution 
time for the given state 


1 + 


/(/ + 1 ) 




( 18 ) 


tR — 


4tt 


-n 




1 + 


1 - 


/(/ + 1 ) 




( 19 ) 


Thus employing as a formation criterion that the bound state constituents life-time r = 1/F^ 
be larger than the revolution time leads us (see Eq. 10) to the inequality: 


(n, 1) Bound state formation 


IOttFo^ 


> 


47r 




n 


1 + 


1 - 


/(/ + 1 ) 




( 20 ) 



We obtain two different conditions for the existence of the IS* and 2P bound states: 

1 S'Bound state formation > 9.15mi, (21) 

and 

2P Bound state formation > 14. 7m^-. (22) 

We thus realize that both formation criteria, Eq. (16) and Eq. (21,22) give comparable 
restrictions on the value of the fundamental scale of SUSY breaking Fq. In order to make sure 
that the bound state can be produced adopting both formation criteria the more conservative 
criterion of the two (that of the revolution time) is chosen. 

For a slepton mass of the order of 100 — 200 GeV the formation of the bound state(s) 
is assured on account of Eq. 7 (v^Eq > A) and the fact that A ^ 10 — 100 TeV. We stress 
that the slepton mass is independent of the energy scale which determines only the 
strength of the gravitino interactions. 

IV. DECAY WIDTH AND ANNIHILATION MODES 

The scalar bound state formed by a pair of slepton NLSP in an e+e“ collision will be a 2P 
state with several decay channels. It will decay into pair of standard model fermions. On the 
other end being formed by NLSP it will have only one decay channel into super-symmetric 
particles: it will decay (annihilate) into a pair of LSP, the almost massless gravitino (gold- 
stino). Finally the 2P state will decay via dipole interactions to a IS state emitting a 
photon. Before discussing these decay channels in detail we make an important remark. 
Within the slepton co-NLPS scenario the bound state can either be ifiifi or rfr)) as in and 
Ti are nearly degenerate in mass. When discussing the decay of the bound state however the 
case of the Ti is somewhat complicated by the fact that left-right mixing must be taken into 
account and diagrams, which are absent in the only case, have to be included. Therefore 
as a first step we consider only a bound state of Ir = cr, JIr. The case of the Ti shall be 
treated on a separate work. 

In addition while p and r cross sections are universal e pair production suffers from destruc¬ 
tive interference effects between s-channel and t-channel diagrams in the parameter region 
relevant for GMSB models (moderate values of tan/?) [14]. For this reason in the following 
we shall consider only the the bound state of JIr (smuonium). 
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The calculation of the partial widths of the decay of the sleptonium bound state 
i.e. 2P —X or IS* —X, is done by relating the amplitude Ai{B X) to that of the 
process X via a non relativistic model of the bound state [15]: 

_ p ^ 1 r / \ 

M(B ^ A-) = yss y — [x {i-A (23) 

iplk) being the Fourier transform of the hydrogen-like wave function of the non relativistic 
bound state. The amplitude A4 xj will in general be described by one or more 

tree-level Feynman diagrams an thus its threshold behaviour s be inferred. 

Thus for each decay one writes down the amplitudes of the contributing Feynman diagrams 
and then extracts the dependence of the full amplitude on k, the momentum of the con¬ 
stituents sleptons which is assumed to be smaller with respect to the mass of the sleptons, 
|fc| -C giving y/s = ^ 

In S'-wave decays one obtains an amplitude whose first term in the momentum expansion 
is a constant. Then the k integration in Eq. (23) gives the wavefunction evaluated at the 
origin: J d^k/{2T[Y'ip*{k) = '^*(0). In P-wave decays one obtains instead an amplitude 
whose hrst term in the k momentum expansion is linear. Then the k integration in Eq. (23) 
gives the gradient of wavefunction evaluated at the origin: J d^k/{27i)^ki r{k) = v^rio). 

For further details see for example [16, 17]. In passing we note that in the process of 
comparing the amplitude of the two photon decay mode of the IS state we hnd an extra 
factor of 2 in Eqs. (Al) and (A2) of ref. [17], relative to stoponium decay into two photons 
and two gluons, a remark that had already been made in ref. [18]. 


A. Decay channels of the 2P state 


Decay into gravitinos: 2P GG 

The process is described by a t-channel exchange of a lepton of the same flavour of the bound 
state. The decay width is: 


F(2P ^ GG) 


1 |i?'2p(0)P 

327r2 M% 


(MA 


(24) 


Annihilation into neutrinos: 2P (£ = e^fi^r). 

This is the simplest decay annihilation channel which takes place only through the Z boson. 


In this case there are no t-channel exchange graphs, even when the flavour of the neutrinos 
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is the same of the sleptonium bound state, since only left s-leptons (f'i) do couple charginos 
and neutrinos. The decay width is: 


r(2P ^ f(r,, e,) (25) 

where f{x,y) = 1/[(1 - + {xyf] and = Mz/Mb] ez = Tz/Mb- 

Annihilation into charged standard model fermions 2P —// 

First we consider the case that / is either a quark or a lepton {£' ^ i), i being the flavour 
of the sleptonium bound state. When this is the case the decay is through the annihilation 
into 7 and Z-boson only: there are no t channel exchange diagrams. The decay width is 
then: 


T(2P^f!) = SCfO- 




m. 

[ —- 

Ml 


mr 


X 


Q} + (1 - 4 ) f{rz, ez) 


A d n J V ■^1 J ^ n • Zj ! J y \ (^6) 

4 COS^ COS^ dYC \ 

where: Cp = ?> for quarks while Cp = 1 for leptons; Q/ is the charge of the fermion in units 
of +e; cy and ca are respectively the vector and axial coupling of the fermion to the Z-boson: 
Cy = Tj — 2QfS,ivfdw and c;^ = T|; m/ is the mass of the fermion; the function f{x,y) is 
the same that appears in Eq.(25). One might notice that the above formula reduces to the 
decay width into neutrinos by taking the limit Qf = 0 and mj = 0. It also reduces to the 
formula of refs.[19, 20]: i?2p(0) being the derivative of the radial part of the wave function 
at the origin, and Mb = 2m^ is the bound state mass. 

When the fermion / is the lepton i of the same flavour of the slepton ip forming the bound 
state Eq. (26) is to be replaced by the following (the lepton i is assumed mass-less): 


r( 2 P^£+r) = 


|fi 2 p( 0 )l 


2 r 


1/2 + 4 sip Sw - 2 sip 6 * 1 + . 

^ + 


^l/2 -2sin^ dw 


{l-rl) f{rz,ez) 


COS^ d\y 

-z -^ OwG — ez) 

2 cos^ dw 


(27) 


where the form factor G describes the diagram of the f—channel exchange of a virtual 
neutralino. Assuming that the neutralino is mostly bino one has the simplified expression 
with: 


G = 


AMl 


Ml + 4m2 
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Dipole decay into the ground state and a photon, 2P —IS* + 7 

The decay to the ground state takes place through a transition with emission of a photon. 
This transition can be computed in the long wavelength approximation. Indeed the photon 
momentnm Q = AE 2 P- 1 S, see Eq. (15), and the bonnd state dimension (Bohr radius) 
rp = 2/{mia) satisfy the relation 

6 

Qrp = —a ^ 2 X 10 ^ << 1 
16 


which is the condition that makes the dipole approximation suitable. Then a standard 
qnantnm mechanics calculation gives [ 21 ]: 


r( 2 P ^lS + ^) = - a{AE 2 -if{D 2 ,if 

y 

where AP 2-1 is the energy of the emitted photon, and 

poo 

D2.1 = ( 2 P|r|lS) = / <irr='Pis(r)P2p(r) 


(28) 


(29) 


is the dipole moment (see [13] and references therein). The wave fnnctions to be used are 
the one of the Coulombic model, given by 




3/2 


Risir) = 2 — exp - 


VbJ 


RipW) — 


(±:\ 


3/2 


rp 


exp 


2rB 


(30) 

(31) 


Kirs) rp 

where rp is the Bohr radius dehned as rp = 2/{mia) = 4/(MbQ;). Using the above wave 
fnnctions one obtains the following expressions for the dipole decay mode: 

64 


r( 2 P ^ 1 ^ + 7 ) = 


p 


6561 


10-2 


(32) 


Then we have to compute the total decay width to fermions: 


r(2P ^ fermions) = ^ r(2P ^ //) + r(2P ^ GG) (33) 

f=qM' 

The derivative in the origin of the radial wavefnnction of the 2P state is easily compnted and 
it follows that: |P 2 p( 0 )P/^b = q;^Mb/ 24576. Inserting the results of the Eqs. (24,25,26,27) 
into Eq. (33) one obtains: 


r(2P ^ fermions) = a^Mp klO ^ 


(34) 
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where k is a numerical constant of order unity with a mild dependence on the mass of the 
bound state Mb, the fundamental scale of SUSY breaking and the mass of the neutralino 
m^o. It then follows that the branching ratio of the dipole decay is Br(2P IS'+y) = 100% 
to within one part in 10^. For all practical purposes the 2P state will decay with probability 
1 to the ground state IS* emitting one photon with an energy of a few MeV. 


B. Decay channels of the IS ground state 


Decay into two photons. 

In this case since the photon is described by transversely polarized states, in the non- 
relativistic limit the t- and u-channel diagrams do not contribute and only the sea-gull 
diagram survives. The decay width is given by: 

.2|fts(0)P 


r(lS ^ 0 - 7 ) = 2 a- 


A/| 


(35) 


Decay into yZ. 

Again the fact that the photon does not have longitudinal polarization states selects, in the 
non relativistic limit, only the sea-gull diagram. The decay width is given by: 

2 siii^Ow |^is(0)P 




Ml 


Ml 

4 


(36) 


cos^ 6w 

Decay into ZZ. 

Here there are four diagrams that give non zero contribution. The diagram with a in in the 
t-channel and its exchange contribute only for the longitudinal polarization states of the Z 
gauge boson. The sea-gull term {iRiR — ZZ) and the s-channel higgs exchange give non 
zero contribution for all type of the Z boson polarization. The decay width is found to be: 

■2 |i?is(o)p. r_ Ml 


r(i^ 


a 


ZZ) = ^ tan'^^w 


-7 


Ml 

I Ml 


(3z{\F + 2\M?>- 


Ml 




4M| 


4M| 


- 1 -2G (3fJeF + 2) 


3M? 


AMI 


8 M| 


(37) 
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where: 


(^z 

F 

G 


Cl 


C 2 


1-4 


Ml 

Ml 


E 

i=l,2 


Ml 


Ml 


Mk 


F ir fjo Mfjo 


4:Ml/{Ml - 2M|) 

cos(/9 — a) cos(/9 + a)/ sin^ 6w 

— sin(/9 — a) sin(/9 + «)/ sin^ 9\y 


while F is a form factor arising from the diagrams with s-channel Higgs exchange, G arises 
from the t-channel £ij-exchange diagrams. 

Decay into W^W~. 

In this case the t-channel exchange of a sneutrino is absent (no z/^j) as well as the sea-gull 
term [Ir Ir — WW), and only the s-channel Higgs contribution is present. The partial decay 
width is: 

Ml I M%' 

where: 



a 


r(l^ ^ H"+H"-) = y 


|flis(0)|- 

Ml 


FfPw 


Pw 


F' 


Cl 


C 2 


T -4 


Mw 

Ml 


Y' ^ _ Mi _ 

cos(/5 — a) cos(/5 -f- a) 

— sin(/5 — a) sin(/3 -f a) 


Decay into hh. 

Within the minimal supersymmetric version of the higgs sector there are hve higgs states: 
three neutrals, (h, H, A) and two charged, {H^). Here we consider only the decay of the IS* 
state into a pair of the lightest Higgs states {h). The process receives contribution from three 
diagrams: a) t—channel £ij-exchange; b) s-channel higgs exchange (h, H) [CP invariance 
forbids s-channel exchange of H]; c) sea-gull term ^R^R — hh. There are no diagrams with 
s-channel exchange of Z since Bose symmetry forbids ZH^H^ couplings. The decay width 
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is found to be: 


r(i^ ^ hh) 


\RlSm% \y 

2 cos Mq ^ 


(39) 


f^h 

Y 

C2 
Cl 

In Fig. 2 the partial widths of the various decay channels of the IS* state are shown with 
respect to the bound state mass Mb- We see that the two photon channel always dominates. 
The branching ratio Br{lS 77 ) is Br = 0.65 at = 100 GeV and Br ^ 0.52 up to 
= 1000 GeV, as shown in Fig. 3. 

The detection of the P wave bound state is therefore associated to the emission of a soft 
photon plus a subsequent emission of two hard photons given by the decay of the scalar 
ground state. 


= ./ 1 - 4 M 
A/| 




M| 


cos(2q;) .0^ .2/ 4M| 

-- —- + sm^ 9w sm^(Q; + /3) 


,tl"2 MI-MI,+iThoMho 2 
= (3/2) sin^(/3 + a) cos(2q;) 

= 2 sin(2Q;) sin(Q; + (3) cos(q; + (3) — cos(2q;) cos^(q; + (3) 


2Ml - Ml 


V. CROSS SECTION 


A. Slepton pair production cross-section 


We use the notation = (3tan6*vK — cot 6w)l‘3 and = (tan 6 'iy + cot6'iy)/4. The 
production cross section for scalar sleptons - except for selectrons - is given by [ 22 ] 


nee 




ttq;" 

T 2 


■/S' 


4 («| + - 4aA', “ ^z) 

~s^ * (s - M|)2 + {MzTzY 


(40) 


with i = L, R] = 2(a^ — (3i) or = 2{ai + Pi) for left- and right- handed sleptons 
respectively. 

It is useful to write down this expression of the P wave cross section in the following 
manner: 

4q;2 ^ 

a{e^e —>• ii) = - P^ d{s) (41) 

3s 


15 



B. Bound state production cross-section 


We shall write the threshold cross section of the bound state in terms of its Schrodinger 
Green function. To review briefly the method, explained in more detail in [23], consider 
the bound state described by a Schrodinger equation with a suitable potential V{x). The 
threshold cross section is then proportional to the imaginary part of derivative taken at the 
origin of the P wave Green function of the problem Gi{x, y,E). E is the energy displacement 
from threshold, and the hnite width of the state is taken into account by the substitution 
E^ E + iT. 

The cross section for the production of a bound state can be normalized to the QED 
process 


a e^e 


aie^e 


Ajr F) f) 

= —d{s)^m Jr —— G^^\x,y,E) 
rrr~ [ oxi oyj 


I 33=0,y=0 


(42) 

(43) 


where d{s) is the usual expression of the Born cross section for the process e+e“ — jly, 
written in Eq. (40). In our investigation the interaction among the two superpartners is 
driven by a Goulombic interaction with 


V(r) = -- (44) 

r 

Adapting the notations of [24] we set E = ^/s — 2m^ as the energy displacement from 
threshold, = —rrij^E, the wavelength A = 9ma/8, and the wave number u = X/k. Now 
the explicit form for the i = 1 Green function Gi of Eq. (43) takes the form 


G,(0,0 ,*:)-|tA{2(i^-V) ± + +27 e-H+v.,( 1-0 + y} (45) 

7 £; is Euler’s constant, is the digamma function, logarithmic derivative of the T function 
and p/ is a soft scale estimated from a relativistic framework (see [23] and references therein). 
The derivative of Eq. (45) needed for computing the cross section Eq. (43) has a simple 
expression, as we have 

Tr ^ ^ G(') {x, y, k) = 9G^^^ (0, 0, k) . (46) 

oXi oyj 3,^0 
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From Eq. (45) one can readily notice that the leading term in the cross section Eq. (43) for 
large E is given by k^, whereas the peaks of the bound state energy levels are determined 
by the digamma function in u. 

The Green function method is a non-relativistic procedure. We have therefore to ensure 
that the velocity of the sleptonium constituents is low enough in order to keep relativistic 
corrections negligible. Starting from the parametrization of the center of mass energy ^/s = 
2mi + E and assuming an upper value for the constituent velocity Pm ax we obtain an upper 
bound for the energy offset from threshold 


E < Em AX — 'i^iPmax 


(47) 


by means of a series expansion in E. This relation translates to the maximal allowed value 
of Lorentz boost parameter 7 

(48) 


Em AX _ 7 MAX 1 


m-n 


'yliAx 


We thus dehne the non-relativistic domain by imposing that the value of 7 differs from 1 by 
less than 10%. This then gives Pm ax ~ 1/v^ and therefore, for a slepton mass mi = 100 
GeV, an energy offset from threshold equal to Em ax = 20 GeV. From eq. 47 one could also 
observe that the acceptable (non-relativistic) threshold energy range increases when larger 
values of the constituent mass are considered. 

In addition we emphasize that the entire procedure of treating the bound state breaks 
down away from threshold, independently of where relativistic corrections may become im¬ 
portant. 

VI. RESULTS AND DISCUSSION 


We shall analyse the threshold behaviour of cross section for a range of masses and SUSY 
parameters for which the bound state formation is envisaged, as discussed in section IV. 
Following [23] we observe that the relevant region for our analysis is the one above threshold, 
i.e. for E > 0. In fact the region below threshold, E < 0, is characterised by peaks in the 
cross section located at the discrete energy values of bound states. Their width is given 
by the annihilation modes which, as shown in IV is of the order of the eV at most. From 
eq. (14) one can estimate the separation of the discrete peaks. They merge when the peaks 
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are not distant enough 


m 


-a 




{n + 1 )^ 


(49) 


The last resolved peak has a quantum number n given by 

2^ + 1 4r ^ 

2 ( 2 I -2 ^ ^ 

Due to the beam energy spread of the collider, much larger than a few eV of the natural 
width of the state and of the order of few GeV [25] this structure cannot be observed. The 
difference from the usual Born cross section results therefore should be sought for i? > 0. 

The effect of the bound state to the cross section at threshold is to accumulate and 
merge the peaks towards the E = Q value, giving a larger result than the naive Born cross 
section, as we could see in Figs. 4 and 5. 


A. The signal and its statistical significance 

Let us now discuss in some detail the signal 

e’'“e“ — 2P —IS* + “soft 7 ” — >• 77 (51) 

where the soft photon is assumed to be undetected since its energy is of only 1 or 2 MeV, 
see Eq. 15. Indeed it is known [26] that the photon energy resolution of the calorimetric 
detector is, for low-energy photons, of the type 5E/E ^ 0.11/ E/GeV which implies that 
5E/E ^ 1 for E = 10 MeV. Therefore the 1-2 MeV soft photon of our signal will surely be 
much below the energy resolution and will be undetected. 

The observed final state is therefore two hard back-to-back photons. This two photon 
signal of the production of the 2P sleptonium bound state is to be compared with the QED 
two photon process 77 which is expected to be the dominant background. The 

number of events of the signal can thus be estimated by: 

Nsign. = klas Br{2P ^ IS +'y) Br{lS 77 ) (52) 

where as is the threshold cross section for the production of the bound state computed with 
the Green function method. Given the fact that Br{2P —IS* -|- 7 ) ~ 1 we have 

Nsign. = Bas Br{lS 77 ) (53) 
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The statistical significance S'S' of the signal can thus be estimated by: 


SS = 


^signal 

^background 


Cas Br{lS 77) 


A comment is in order at this point. The two-photon decay of the bound state is isotropic 
in the rest frame of the bound state which however will approximatively coincide with 
the laboratory frame as we are assuming threshold production with E < 20 GeV. The 
distortion introduced in the laboratory system when boosting the isotropic distribution of 
the two-photon signal from the bound state rest frame depends on the j3{M) of the bound 
state. Indeed in the bound state rest frame let 6* be the angle of the direction of the two 
outgoing (back-to-back) photons: 

d(TB _ CTsis) . X 

dcose* ~ 4 ^ ’ 

where the phase space integration has been reduced of a factor of 1/2 due to the identical 
particles in the final state. The distribution can be boosted [27] to the laboratory frame 
where the bound state has a velocity P{s), ^/s = 2m^ + E, E being the offset from the 
threshold: 

d(TB _ (^b{s) 1-/3^ 

dcos6 4 {l — /3cos6y 

where now 6 is the direction of the photon in the laboratory frame. 


On the other end the angular differential distribution of the photons for the QED process 


77 is given by: 


dcos9 


ttq;^ 1 -|- cos^ 6 
s 1 — cos^ 6 


which is symmetric and peaked in the forward and backward directions. Again a factor of 
1/2 has been included in the phase-space factor due to the identical particles of the final 
state. Applying an angular cut in the forward and backward directions 

Of) < 6 < TT — Of) 


it is possible to suppress the background in such a way as to obtain an interesting statistical 
significance. We define a cut dependent statistical significance: 


^^(^ 0 ) = 


NsigiOo) _ Br{lS ^ sin6'd6' {daB/dcos9) 


'N,ack{9o) 


'C jg sin 9d9 [do^f^jd cos 9) 
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The statistical significance as function of the angular cut is shown in Figs. 6 and 7. We 
observe a very mild dependence on 9q apart from values of 6^0 ~ 0 and do ~ 7 r/ 2 . For these 
limiting values the statistical significance approaches zero. When 9o —^ 0, SS{9o) —0 since 
the QED two photon background total cross-section is divergent. On the other hand when 
9o ^ 71/2 one is integrating over a region of phase space of vanishing measure so that both 
cross sections vanish as Zq = cos^o in this limit. One notices that for values of 9o in the 
range 0.4 < 9o < 0.8 the statistical significance is almost constant and ranges from SS ~ 4 
to 5'S' ~ 11 for values of the energy offset respectively of = 10 GeV and E = 20 GeV, 
and for a fixed slepton mass = 100 GeV. 

In Fig. 8 we show the dependence of the statistical significance SS{9o) computed with 
9o = 0.55 which corresponds to the maximum of the curves in Fig. 6 for two given values of 
the energy offset E = 10 GeV and E = 20 GeV as function of the sleptonium bound state 
Mb = We see that even going to higher values of the slepton mass ~ 200 GeV 

the statistical significance is still at an interesting value SS ^ 2 (at = 20 GeV). 

It should also be checked that along with interesting values of the statistical significance 
one has also interesting absolute values of the cross sections, i.e. an observable number of 
events. Indeed from Fig. 4 at = 20 GeV, = 100 GeV, one can see that as x BR{ 1 S —> 
77) ^ 80 X 0.63 fb ~ 50 fb. This would correspond to a total number of signal events 
Nsig ~ 5000 assuming an annual integrated luminosity Lq = 100 fb“^. For E = 20 GeV, 
= 200 GeV one finds as x BR{ 1 S 77) 8 x 0.63 fb ~ 5 fb. This would correspond 

to a total number of signal events Ngig ~ 500, ten times lower than the 100 GeV mass case. 
The above estimates have been done using the total cross section as- We have checked that 
taking into account the distortion of the boost to the laboratory system from the rest frame 
of the bound state, where the two-photons are isotropically distributed, changes the above 
estimates by an amount which is at most of 1 %. 
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VII. CONCLUSIONS 


In this work a novel signatnre at the next LC is proposed for the detection of a sleptonium 
bonnd state (of a charged NLSP) assuming that the energy of the collider happens to be 
around the threshold. It is well known that GMSB models are characterized by large regions 
of the parameter space in which the NLSP is a charged slepton = (cr, ri). While 
interesting signatures of the production of a pair of charged sleptons NLSP have already 
been considered in the literature, here we discuss the formation, decay spectroscopy and 
possible detection of the NLSP (charged slepton) bound state and its possible signature at 
the next LC. For the sake of simplicity we have considered only the case of = JIr. 

At an collider the smuonium bound state is produced in a 2P state. We assume 
a standard Coulombic interaction to describe the bound states and thereby estimate the 
energy levels and their separation which are needed to establish whether or not the formation 
takes place. As a criterion for formation we have chosen the more conservative condition, 
as opposed to the level gap criterion, obtained by requiring that the revolution time be 
larger than the lifetime of the rotating constituents. We find out that the formation of the 
2P bound state is assured for > 1A.7mi and thus for all relevant values of V^, the 
fundamental scale of SUSY breaking. We have first analyzed the decay channels of the 2P 
state. These include: a) decay into a pair of gravitinos, 2P —> GG, b) decay to a pair of 
standard model fermions 2P ff and c) the dipole transition 2P IS* + 7. Studying the 
partial widths and the branching ratios we find that for all practical purposes the 2P state 
decays to the IS* state emitting a photon whose energy is of a few MeV (AP 2 P-is)- 

In turn the decay channels of the IS* state are to the following final states: 
'y'y,'yZ, ZZ,W^W~ ,hh and it turns out that the dominant decay is IS —> 77 as shown 
in Figs. 2,3. Thus as the MeV photon from the 2P decay goes undetected because its energy 
is below the detector resolution, the following signal of the sleptonium bound state is defined 

e’'"e“ — 2P —IS* + “soft 7” —77 , 

where the observed final state is two hard hack to back photons, and ’’practically” no missing 
energy. 

The bound state production cross section has been estimated using the Green function 
method. Due to the fact that there a bound state actually exists, this effect accounts for 
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a dramatically different cross section with respect to the case of the Born production. For 
energies just under the threshold, there are several peaks centred at the discrete energy 
levels of the bound state. Approaching the threshold level from below, those peaks merge 
and accumulate towards E = Q level. This fact translates to the energies above threshold 
as well {E > 0), increasing the cross section of the continuum. The net effect is a larger 
cross section with respect to the naive Born case. This effect grows with the strength of the 
coupling of the particles which form the bound state itself. It is still very noticeable even 
for a weak coupling like the sleptonium case, as it is larger than the Born cross section for 
about 10% to 30%, the difference increasing with the displacement from threshold. 

We have given an analysis of the two-photon signal comparing it to the QED two-photon 
cross-section (at leading order) dehning a statistical signihcance depending on the angular 
cut which is introduced in order to reduce the QED background. We hnd that at an energy 
offset of E = 20 GeV from the threshold, the statistical significance SS{6q) (computed 
for an angular cut 9q = 0.55 which maximizes it) goes from SS = 11 at = 100 with 
Nsig ~ 5000 to SS = 2 for = 200 with Wig ~ 500. 

Our study has been done assuming a charged slepton co-NLSP scenario which is a pecu¬ 
liarity of GMSB models. On the basis of the results discussed above we conclude that the 
study of a charged slepton NLSP bound state through the two-photon signature at the next 
LG has the potential to shed insights into the mechanism of SUSY breaking. 

One of the authors (N.F.) wishes do dedicate this work to the memory of Edmondo Pedretti. 
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Slepton co-NLSP scenario 




E 



FIG. 1: Sparticle spectrum in the slepton co-NLSP scenario. Only the masses of the sparticles 
which are relevant for our processes are shown as a function of the NLSP mass rriNiSP = 
What is shown here is a scan of the GMSB parameter space done imposing the slepton co-NLSP 
scenario condition given in Eq.(9). 
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FIG. 2: Partial widths of the various decay channels of the 15 state, as a function of the bound 
state mass Mb = 2mMLSP- 
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FIG. 5; Comparison of the Green function method and the Born level expression for the cross 
section, for the right smuon. Here we assume that = 200 GeV and a displacement from 

threshold of up to 20 GeV. 
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FIG. 6 ; Statistical significance (signal to noise ratio) of the cross section given by the Green 
function method and the 77 QED background as a function of the angular cut angle 9o, for the 
right smuon /i/j. We assume = 100 GeV; the various curves represent different displacements 
from threshold. 
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FIG. 7: Statistical significance (signal to noise ratio) of the cross section given by the Green 
function method and the 77 QED background as a function of the angular cut angle 9o, for the 
right smuon /i/j. We assume = 200 GeV; the various curves represent different displacements 
from threshold. 


32 









FIG. 8; Statistical significance (signal to noise ratio) with respect to the mass of the bound state 
Mb = 2mNLSP- We show here the case of the right smuon fiR bound state produced at an energy 
£' = 10 and i? = 20 GeV above threshold. The angular cut is assumed to be at 9q = 0.55, maximal 
value of the signal to noise ratio. 
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